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Abstract 

A numerical investigation of the quenched Schwinger model on the lattice 
using the overlap Dirac operator points to a divergent chiral condensate. 
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I. INTRODUCTION 



Use of quenched QCD as an approximation to the full theory depends upon a good 
understanding of the regions of parameter space where the quenched theory differs in im- 
portant ways from the full theory. For the case of the chiral condensate (ipip) , there may be 
qualitatively different behaviors for sufficiently small quark mass. Whereas the condensate 
is expected to be finite in the full theory, there are theoretical arguments |l| and some initial 
numerical indications B that it diverges in quenched QCD. Also numerical analysis || of 
an instanton gas model shows a divergence. A careful study, using a lattice Dirac operator 
that obeys chiral symmetry, to determine the mass range in which quenched QCD is a good 
approximation to the full theory has not been done. 

Similarly for the Schwinger model, the full theory has a finite condensate M but there 
are predictions |5],|6j that it diverges in the quenched theory. Thus the Schwinger model can 
be used to investigate the two-dimensional versions of these questions concerning anomalies, 
topology, chiral symmetry, and condensates and their impact on the relationship between 
full and quenched theories. Although the Schwinger model is in most ways much simpler 
than QCD, it does present some peculiar difficulties of its own. Strong infrared effects, which 
are dynamical and nonperturbative in QCD, are already kinematical in the lower dimension 
of the Schwinger model. There is the possibility that infrared enhancement in the quenched 
case gives a fermion spectral density that is divergent as the eigenvalue A goes to zero and 
that there is a corresponding infinite condensate < ipip >■ We have investigated this issue 
numerically using the overlap Dirac operator to describe the massless limit for the fermions 
and have found strong evidence for these divergences in the quenched Schwinger model. 

When stated in terms of the low eigenvalue behavior of the fermion spectral density in 
the quenched Schwinger model, theoretical discussions have given a lower bound that is 
finite |7[ and stronger one that is divergent 0. Some estimates |^J,|6| that are not bounds 
have suggested a form diverging exponentially in the volume e° 9 . The discussions in 
Refs. PJ^] are given in terms of the eigenvalue shifts of the would-be- zero modes associated 
with subregions of the whole two-dimensional volume V. With larger shifts due to 
interactions with other subareas, the spectrum is flat at small A in the infinite volume 
limit. Smaller shifts leave the would-be- zero modes concentrated near the origin so that 
the spectral density there diverges as V — > oo. Other arguments HQ for the form of the 
divergence proceed along different lines, but the implication for the spectrum is that the 
lowest eigenvalues are exponentially small in the volume with a corresponding exponentially 
large spectral density and condensate. 

The data that we present here covers a range of lattice sizes from 8 2 to 32 2 . The full 
spectrum of the overlap Dirac operator was calculated in gauge backgrounds from the Wilson 
action at several bare couplings. The next section discusses the lattice formalism used in 
this paper. The third section discusses the physics issues in more detail. The fourth section 
gives our numerical results. In addition to the spectrum itself, there are measures of its 
behavior including < ipip > and the distribution of the lowest eigenvalue. The last section 
contains a summary of our results and some concluding discussion. 
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II. LATTICE FORMALISM 



Since we are interested in studying the small mass region and the massless limit, we need 
to work with a lattice Dirac operator that respects chiral symmetry. We will use the overlap 
Dirac operator for our numerical study. It has the form [H] 



2 



1 + m + (1 - m) lh e{H w ) (1) 



with H w being the hermitian Wilson Dirac operator in the supercritical region and < m < 1 
is the bare fermion mass. The hermitian overlap Dirac operator H = 75D has paired non-zero 
eigenvalues. The topological zero modes are chiral and have partners with unit eigenvalue 
and opposite chirality. In a fixed gauge field background 0, 

j = \Q\ 1 v 2m(l - A 2 ) 
<W> mV V^ \ 2 (l-m 2 )+m 2 ' U 

The sum is over all positive non-zero eigenvalues of H, Q is the global topological charge, 
and V is the lattice volume. 

In the numerical calculation, we generate gauge fields distributed according to the Wilson 
gauge action 

s* = i E RqU p ( 3 ) 

9 p 

with Up the product of U(l) link elements around a fundamental plaquette and g the lattice 
coupling constant. The fermions have periodic boundary conditions []. For each choice of 
g and L, we diagonalize H w in a fixed gauge field background and form H by first forming 
e(H w ). We then diagonalize H 2 in the chiral sector that contains topological zero modes, if 
any. Since all computations are done in double precision, we know the non-zero eigenvalues 
of H to an absolute precision of 10~ 8 . In addition we know the exact number of zero 
eigenvalues of H by counting the difference between the number of positive and negative 
eigenvalues of H w [jTH . 



III. PHYSICS ISSUES 

In the multi-flavor Schwinger model, the classical U(l) chiral symmetry is explicitly bro- 
ken by the anomaly, while the SU(N) chiral symmetry cannot be broken in two dimensions. 



This choice of boundary conditions is not as restrictive as it seems since we only have one 
fermion. A gauge field configuration can be multiplied by an arbitrary constant U(l) field on each 
link in either of the directions without changing the gauge action. Since all of these possibilities are 
included in the sum over gauge field configurations, there is no real distinction between periodic 
and anti-periodic boundary conditions. More generally all boundary conditions that are periodic 
up to a phase are equivalent. 
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The 't Hooft vertex < fli V'iV'i > is not associated with an intact symmetry or Goldstone 
bosons, so it can and does have a nonzero value JO]]. In the quenched case, the exact zero 
modes of the massless Dirac operator cause a divergence in < ipip > in the massless limit 
at finite volume. But as seen in (§), the divergence is of the form < \Q\ > /(mV). Since 
< \Q\ >oc W, it follows that 

lim lim < \Q\ > l(mV) = 0. (4) 

m— >0 V— >oo 

This trivial divergence does not contribute in the case where one first takes the thermo- 
dynamic limit and then takes the massless limit. But this finite volume divergence does 
not appear in the unquenched Schwinger model. The zero modes of the Dirac operator in 
these backgrounds cause a suppression of such gauge field configurations when the fermion 
determinant is included as part of the gauge field measure. 

The small eigenvalue behavior of the spectrum determines the contribution that the 
second term in (|2[) makes to < ipip >■ Thus the issues to be numerically investigated are 
centered upon the small eigenvalue behavior of the massless Dirac operator. The main 
question is whether the infinite volume spectrum is flat as the eigenvalue A goes to zero or 
has a divergence at small A. 

Consider the gauge field seen by the fermion. The plaquette magnetic field of the 
quenched Schwinger model is ultra-local with the field fluctuations on different plaquettes 
uncorrelated in infinite volume. The plaquette angles are approximately gaussian distributed 
at weak coupling. Thus the study of fermionic observables in the quenched theory is best 
thought of as an investigation of a disordered system |7|]. 

Let us begin the discussion with two much simpler examples. For the case of free fermions 
on an L x L lattice with periodic boundary conditions, the low-lying levels are 

A^[(^i) 2 + (^) 2 ] 1/2 (5) 

so that the level spacing is of order 1/L, and the density of states per unit volume is of order 
A at the low end. 

Another simple case is a uniform magnetic field B, which gives Landau levels. The level 
spacing is of order B, and the degeneracy of each level is of order BV . With the scale of 
energy intervals larger than B, the density of states is flat. As we will see later, the typical 
BV is gL so that the average of B over V does get smaller with increasing volume. Thus 
the Landau levels give a flat spectral density if the energy resolution is coarser than g/L. 

For the case at hand of particles with gyromagnetic ratio 2, there is a cancellation between 
the paramagnetic magnetic moment interaction with the field and the diamagnetic kinetic 
energy contribution that puts the lowest Landau level at exactly zero energy. These are the 
BV I (2ir) zero modes. 

When the net flux is zero and all boundary conditions are periodic so that the vector 
potential can be put in the form 

A* = e n^(f>, (6) 
there is also a pair of zero modes with opposite chirality. These have the form 
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iP + = e' 




1 




) 



and ip. 



'_ = e 




(7) 



For the quenched Schwinger model the field is neither zero nor uniform. As noted above, 
the magnetic field is random with no plaquette-plaquette correlation between the magnetic 
field on different plaquettes. Thus the variance increases as the area. With the coupling g 
small and gR large, the flux through the area R x R is gaussian distributed with a typical 
size of gR, so that the area average of the field strength is B = g/R. 

What is the fermion spectrum in that case? The index theorem tells us that for a net 
flux 27r/ through the area L 2 , there are / modes with zero eigenvalue. For / = and all 
boundary conditions periodic, there are two zero modes of the form above. But what else 
happens at the low end of the spectrum? There are two suggestions for an answer in the 
literature. The discussion of Casher and Neuberger J7j begins by dividing the L x L volume 
into R x i?-sized pieces with gR and L/R large. Considered in isolation, each of these areas 
has of order gR zero modes. It is then argued that the effect of interaction between different 
regions is to shift the eigenvalues of the zero modes away from zero, in such a way as to 
produce a spectrum that is bounded below by one that is flat at small A. 

With a similar approach, Smilga || argues for a stronger lower bound that gives a spectral 
density that diverges as A — > 0. His stronger result follows from using the fact that the value 
of the zero mode wave function on the boundary that separates a region with magnetic field 
from one with zero field is exponentially small in the flux through the region. Arguments 
using other methods in his paper and in the paper by Diirr and Sharpe || conclude that 
the divergence is quite strong with a factor e cg v . This would be a consequence of modes 
with eigenvalues as small as the inverse of that factor. 

For a numerical test of the argument used by Smilga to produce the stronger bound, we 
construct a background gauge field configuration that has two regions of size R 2 each with 
constant magnetic field and opposite net fluxes of magnitude R. We study the low lying 
eigenvalues of the overlap Dirac operator and show that they go down exponentially with 



On the L 2 lattice, fix two regions of size R 2 separated by (L/2 — l,L/2 — 1). The slightly 
off-symmetric separation is chosen to avoid any accidental lattice symmetries. On one region 
of size R 2 , we make up a constant magnetic field of flux R, and on the other region, we make 
up a constant magnetic field of flux —R. Elsewhere the field is zero. An initial numerical 
check with the field set to zero in one of the regions, verified the presence of R topological 
zero modes for the overlap Dirac operator. Then returning to the case of interest with 
the field in both regions, we calculated the spectrum again. In Fig.[I], we plot the low end 
of the positive half of the spectrum as a function of R. The lowest few of these small 
eigenvalues go down exponentially in R. We verified that this behavior remains unchanged 
when small random perturbations are added to the link elements. The numerical results 
are less restrictive than the theoretical arguments in that they do not rely on a variational 
argument, which without further work, applies only to a single mode on the lattice. Also all 
the lattice modes in addition to the would-be-zeros are included. This confirms the crucial 
point in the argument for the stronger bound of Ref. M. However, it does not provide 
evidence that the spacings might be as small as e~ cg v . 

Let us now consider the continuum limit. Given the remark above on the lack of cor- 



R. 
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FIG. 1. Low end of the positive half of the spectrum on a L = 32 lattice on the special 
configuration as a function of R which defines the area of constant magnetic field of flux R. 
Eigenvalues below 10 -7 are not shown in this plot. Only for R = 3 are all the eigenvalues above 
this bound. The lines connect eigenvalues at the same position in the ordering for each R. 




relation in the field strengths, it is not possible to use a scale from that in defining the 
continuum. However, there is another simple approach. For g small, ask how large does R 
have to be so that the typical flux through the R x R area is of order one? Since the flux 
variance for a single plaquette is g 2 and the uncorrelated fluxes add, the variance for the area 
is g 2 R 2 . Thus the typical flux through the area is gR, and g sets an inverse length or energy 
scale. This happens to be of the same order as the scale for the unquenched Schwinger 
model, in which the mass in lattice units is g/ s/k. We may hope to get a sensible continuum 
limit by measuring continuum dimensionful quantities in units of appropriate powers of g/a. 
To get the finite volume continuum limit, we will want to take g to zero and L to infinity 
with gL fixed. Lattice eigenvalues, < ipip >, and other quantities with continuum units of 
energy should be considered in ratios like X/g as g — > 0. 

Finally, let us discuss the range of g and L where these interesting effects might be seen. 
From two points of view, we can see that gL must be large. First if L is fixed and g is 
small, then there are only perturbative effects from the gauge field, and the strong infrared 
fluctuations cannot appear. Also if gL is small, then there are essentially no would-be-zero 
modes that could realize the the physical pictures of [0] and ||. If g is large, then the 
gauge field is very rough on the scale of a single lattice spacing, and the continuum-based 
arguments above do not apply. The smallest region that typically contains a unit of flux 
should be several lattice units so that it is large enough for the fermion to realize zero modes 
from the non-zero flux. Thus we must have g small and gL large, which means, of course, 
that L must be large. 

IV. NUMERICAL RESULTS 

The numerical results described in this section give substantial evidence that the infinite 
volume limit of the spectral density p(X) is indeed infinite for A — > 0. We will show this 
by computing the spectrum of the overlap Dirac operator in U(l) gauge field backgrounds. 
The massless limit is approached in the conventional way by adding a standard mass term. 
We will also compute < ipip > as a function of mass in finite lattice volumes and show that 
it grows as the mass is lowered before it finally dives to zero as it must for m = and finite 
volume. We will also show that the average value of the lowest eigenvalue does not scale 
with the volume, nor does it fit predictions from chiral random matrix theory. 

Strong evidence for a divergence in the non-topological piece of < iptp > is seen by 
plotting the gauge ensemble average of the second term in (0) as a function of m at a fixed 
g = OAy/n for several lattice sizes. We focus on the small mass region in Fig.^. The data 
at L = 32 show a rise in < ipip > as the mass is decreased. This divergence is due to an 
accumulation of very small eigenvalues at larger L as seen in the histogram of the small 
non-zero eigenvalues in Fig|| Even though only L = 32 shows a rise in < ipip > at small 
masses, an anomalous accumulation of very small eigenvalues is evident at L = 24 in Fig|| 
The accumulation is not enough to give a rise in < > on the L = 24 lattice. 

The smallest eigenvalue A m ; n has to scale like 1/V for a finite value of the density of 
eigenvalues at zero p(0) and a finite value of < ipip > in the massless limit. In Fig.^, we 
plot the histogram of A m i n L 2 for the various ensembles in the zero topological sector. We see 
that there is no evidence for scaling in the distribution, whereas chiral random matrix theory 
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FIG. 2. Plot of < ipif> > with respect to m at g = 0.4^/7? for four different lattice sizes 




FIG. 3. Plot of the distribution of the small non-zero eigenvalues at g = OA^/tt for four different 
lattice sizes 
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FIG. 4. Plot of the distribution of the smallest eigenvalue scaled with the volume and restricted 
to the zero topological sector at g = OA^/tt for four different lattice sizes. 
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FIG. 5. Ensemble average of the smallest eigenvalue in the zero and unit topological sectors 
at g = OAy/n as a function of the lattice size and plotted to facilitate comparison with the simple 
functional forms described in the text. 
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predicts a universal function of the form |e~~ with z = £L 2 A m i n and £ the value of the 
chiral condensate. A previous analysis fT2| of the distribution of the low lying eigenvalues 
done at smaller physical volume showed that the distribution did not fit the predictions of 
unitary chiral random matrix theory ||13|| . In Ref. [|12] this was attributed to finite volume 
effects. In our case, the reason for the discrepancy is not small volumes but a divergent chiral 
condensate. Now consider the average of the smallest nonzero eigenvalue as a function of 
L. Three simple functions motivated by heuristic physics arguments are c/V, e~ cL , and 
e~ cV . To determine which of these forms is closest to the data, we have plotted in Fig.|5| 
K(A min )/16, ln((A min )/(A min ) |l= 8 )/(8 - L), and 50 ln((A min )/(A min ) |l= 8 )/(64-K) verses L 
for Q = and \Q\ = 1. (The normalizations are just for convenience.) To the extent that 
one of these functions represents the data well, the corresponding graph in Fig.|5| should be 
flat. Evidently e~ cL is preferred. Recall that this is the form that appears in the argument 
for the lower bound in || and in the spectrum from the artificial configurations described 
in Section III. 

Fig.| shows that L = 32 is needed at g = QA^Jn to see the divergent behavior in the 
chiral condensate. This corresponds to a physical volume of gL/y/n = 12.8. To study 
the effect of lattice spacing, we compared this result with others obtained on L = 24 at 
= 0.4(32/24) and on L = 28 at g/y/n = 0.4(32/28). These have the same physical 
volume gL/y/n = 12.8 but are coarser lattices. The comparison in Fig|| shows that the 
divergence visible on the L = 32 lattice is not seen on the L = 24 lattice. However, the 
L = 28 data follow the L = 32 data to very small masses and into the region where the 
condensate begins to grow. We have used dimensionless quantities in this plot to facilitate a 
proper comparison. The scaling behavior is good until m/g gets small enough to emphasize 
the very smallest eigenvalues, some of which are being distorted on the coarser lattice. At 
smaller coupling and at the same physical volume, the scaling behavior extends to smaller 
m/g. 

To illustrate the point that at fixed L, g can be neither too big nor too small if the small 
A growth is to be seen, we have data from L = 32 and four couplings in Fig.|7]. The smallest 
value of g, which corresponds to a physical size of 9.6, shows no growth at all. The medium 
values at sizes 12.8 and 16 show the effect. (Note that the finite size effects between these 
two are small.) The largest value of g has size 32 but the gauge field there is too rough, and 
the small A peak is gone. 



V. CONCLUSIONS 

We have shown that for small coupling and large volume, a small eigenvalue peak appears 
in the spectral density of the overlap Dirac operator and in < ipip >. This is strong evidence 
for the predictions |||6| that these quantities diverge in the infinite volume limit of the 
quenched Schwinger model. There is some evidence that the divergence could be as strong 
as e cgL , but the lattice sizes are insufficient to provide evidence for the stronger e dg2L2 
predictions. Similarly there is limited evidence that the would-be-zero modes of subregions 
of the lattice can provide a physical understanding of the results. However, a definite test 
of that model also awaits data from larger lattices. 

The results in this paper clearly point the direction for further work. The calculations 
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FIG. 6. Plot of < ipip > /g with respect to m/g at gLj^pK = 12.8 on L = 24,28,32. 




FIG. 7. Plot of < i/ji/j > jg with respect to m/g at L = 32 and four different couplings. 
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should be extended to larger lattices so that there are several sizes showing the small A 
growth of the spectral density and < ipip >. With that, it would be possible to study the 
volume dependence of the small A peaks and test in more detail the theoretical expectations. 

Although the quenched Schwinger model is quite some distance from full four-dimensional 
QCD, results from it help to map out the range of territory available to massless fermions 
responding to a gauge field. 
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